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ABSTRACT 

CoRoT-7 b is the first confirmed rocky exoplanet, but, with an orbital semi-major axis 
of O.Of 72 AU, its origins may be unlike any rocky planet in our solar system. In this 
study, we consider the roles of tidal evolution and evaporative mass loss in CoRoT-7 
b's history, which together have modified the planet's mass and orbit. If CoRoT-7 b 
has always been a rocky body, evaporation may have driven off almost half its original 
mass, but the mass loss may depend sensitively on the extent of tidal decay of its orbit. 
As tides caused CoRoT-7 b's orbit to decay, they brought the planet closer to its host 
star, thereby enhancing the mass loss rate. Such a large mass loss also suggests the 
possibility that CoRoT-7 b began as a gas giant planet and had its original atmosphere 
completely evaporated. In this case, we find that CoRoT-7 b's original mass probably 
didn't exceed 200 Earth masses (about 2/3 of a Jupiter mass). Tides raised on the host 
star by the planet may have significantly reduced the orbital semi-major axis, perhaps 
causing the planet to migrate through mean-motion resonances with the other planet 
in the system, CoRoT-7 c. The coupling between tidal evolution and mass loss may be 
important not only for CoRoT-7 b but also for other close-in exoplanets, and future 
studies of mass loss and orbital evolution may provide insight into the origin and fate 
of close-in planets, both rocky and gaseous. 
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1 INTRODUCTION 

The recently discovered exoplanet CoRoT-7 b has a mass 
of M v = 4.8 ± 0.8 Earth masses (MEarth) and a radius 
of only R p = 1.68 ± 0.09 Earth radii (R Ear th). Thus its 
density likely lies between 0.71 and 1.4 times Earth's den- 
sity, which suggests CoRoT-7 b is made o f similar mate- 
rials (|Leger et al.l |2009i ; iQueloz et al.l 120091 1 . It appears to 
be the first example of a rocky exoplanet. However, with a 
semi-major a xis a of only 0.0172 ± 0.00029 AU and e « 
IjLeger et al . 2009), CoRoT-7 b has an orbit unlike the rocky 
planets in our solar system. CoRoT-7 b's surface may even 
be par tially melted as a result of the strong stellar inso- 
lation (|Leger et al.l 12009). (Consequently, when we refer to 
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CoRoT-7 b's "rocky" surface, we are referring to the com- 
position and do not mean to imply that the surface is solid 
rock.) The announcement of an additional, nearby planet in 
the system, CoR oT-7 c, (M v = 8.4 ± 0.9 M Ea rth, a = 0.046 
AU, and esiO) (JQueloz et al.ll2009f ) also suggests that the 
system has an interesting dynamical history- 
Missions like K epler (JBorucki et al.ll2008h and CoRoT 
( Bagli n et al.ll2006r ) may find dozens of additional rocky ex- 
oplanets through transit observations. Although Kepler is 
expected to find several ro cky planets in orbits like the 
Earth's (|Borucki et al.l 120081 ). the first and most easily de- 
tectable rocky plane ts will have orbits closer to their stars 
(Dom ing et al.l l2009b Thus a detailed understanding of such 
planets will be critical to understanding the planets first dis- 
covered by these space-based missions. 

The origins of these planets may be very different from 
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the origins of rocky planets in our solar system. It may be 
difficult for rocky exoplanets to form as close to their stars 
as CoRoT-7 b and c because high temperatures in the proto- 
planetary gas disk m ay inhibit condensat ion and accretion of 
solid materials fe.q.. |Pollack et al.lll994T ). At the very least, 
solid planets that form close-in will probably be dep leted 
in volatile materials (e.g., water [Raymond et al.ll2008l |V In- 
stead of forming in their current orbits, close-in rocky exo- 
planets like CoRoT-7 b may have coalesced in orbits farther 
away from their host stars, where disk temperatures were 
lower, and were brought to their current orbits. A variety of 
processes for br inging planets to such close-in orbits have 
been proposed JMardling fc Linl |2004| ; iGaidos et ail 120071 ; 
iRavmond et alj 1200a ) . and the orbital architecture of sys- 
tems with close-in planets and the planets' physical prope r- 
ties provide clues to their histories (jRavmond et al.l l2008). 

One process that may have been important during 
CoRoT-7 b's history is evaporative mass loss. Evaporation is 
thought to influence gaseous exoplanets made primarily of H 
and He. For example, observations of HD 209458 b suggest 
the planet's atmosphere is evaporatin g at a rate of order 
0.1 M Earth /Gyr |Vidal-Madiar et all 12003. 2004), in line 
with theo retical expectations from studies of evaporative 
mass loss (Lanimer et al. 200 31: iLecaveli er des Etang s et al] 
|2004| ; lYellell2004l ; iTian et al]|2005l ; lYelle et al.ll2008l ) More- 
over, gaps in the distributions of orbital and physical prop- 
erties of close-in planets have also been interpreted to in- 
dicate that many gaseous plan ets have had their atmo- 
spheres completely evaporated (jDavis fc Wheatlevl l2009n . 
The complete removal of a gas giant's atmosphere would 
likely leave behi nd its rocky core, with a mass o f perhaps 
several M Ear th l|Baraffe et alj|2004l ; [Raymond et~afll2008h . 
These considerations suggest there may be a large popu- 
lation of close-in planet ary bodies that are remnant core s 
of evaporated gas giants (jBaraffe et al.ll2005l ; lHattoril l2008'). 
On the other hand, other studies argue that the observa- 
tions o f HD 209458 b do not imply sign ificant loss of mass 
l|Ben-J affel 2007; Holmst rom et al.ll2"008T ). and a theoretical 
study bv lMurrav-Clav et al] |2009) suggested that complete 
evaporation of a gas giant's atmosphere is unlikely. These 
competing hypotheses may now be testable. With the detec- 
tion capabilities of the Kepler and CoRoT missions, rocky 
planets arising from a variety of histories may be detected. 

Evaporation of mass may continue even after the planet 
loses its original H/He envelope. Valen cia et al] (J2009) sug- 
gest the strong insolation received by CoRoT-7 b may be suf- 
ficient to have evaporated and removed several Earth masses 
of material, i f the planet has always b een solid (made of ice 
and/or rockV lSchaefer fc Feglevl (J2009J) showed that thermal 
evaporation of rocky material from CoRoT-7 b's surface may 
produce a tenuous exosphere, similar to Mercury's. In this 
study, we consider evaporation both of a putative gaseous 
envelope and of solid material during CoRoT-7 b's history, 
including effects of its evolving orbit. 

In addition to mass loss, tides have played an im- 
portant and interrelated role in CoRoT-7 b's history. 
Tides have shaped the distribution of orbit al elements 
for c l ose-in exoplanets, re d ucing eccentricities (|Rasio et al] 
1 19961 ; iFord fc R asio 2006; Jacks on et ail l2008h and semi- 
major axes (Barker fc Ogilviei 120091 ; lJackson et al] 120091 ; 
iLevrard et al] 12009). Tides may also help to br ing rocky 
planets into close-in orbits IjRavmond et al. 2008). In fact, 



going forward in time, CoRoT-7 b may migrate into its host 
star in less than a few billion years, and it may ha ve un- 
dergone past orbital migration II Jackson et al] 120091 ) . Also, 
CoRoT-7 b likely has undergone or currently undergoes 
strong tidal heating, either because its original orbit was 
eccentric or interacti ons with nearby C oRoT-7 c keep the 
eccentricity non-zero IjBarnes et al]|2010l ) . Tidal evolution of 
its orbit may have been even faster in the past if CoRoT-7 b 
once had a massive gaseous envelope because the rate of evo- 
lution of a circular orbit scales with its mass. As the planet 
lost mass, though, the tidal evolution may have slowed as 
the planet reached its current orbit. 

In this paper, we study the coupling of evaporative mass 
loss and tidal evolution on CoRoT-7 b. Using parameter- 
ized models for mass loss and orbital decay due to tides, 
we evolved CoRoT-7 b's mass and orbit backward (and for- 
ward) in time to determine the range of original masses and 
orbits consistent with the planet's current mass and orbit. 
Although the exact initial conditions depend sensitively on 
the chosen model parameters, we considered a wide range 
of parameters to map out all the possibilities. We find that, 
if CoRoT-7 b has always been a solid planet, it may have 
lost as much as half its original mass and had its semi-major 
axis reduced by as much as 60% by tidal migration alone. 
If CoRoT-7 b was originally a gas giant planet and only 
recently lost its gaseous envelope, its original mass did not 
exceed 200 MEarth, and its original semi-major axis could 
have been almost twice as large as the current value. In both 
scenarios, CoRoT-7 b may have encountered and migrated 
through a 3:1 mean motion resonance with nearby CoRoT-7 
c, as tides caused planet b's orbit to decay. Consideration of 
the resulting dynamical interactions may place constraints 
on the orbital history of the CoRoT-7 system and the orig- 
inal mass of CoRoT-7 b. We also discuss how the coupling 
between orbital decay and mass loss may play an important 
in the evolution and survival of gaseous close-in planets. 



2 MODELS 

Several recent studies of evaporation of exoplanetary 
atmospheres provide si mple, parame t erized models for 
mass loss rates (e . g.. [ Baraffe et al] |2004| ; Ijaritz et al] 
120051 ; lErkaev et all l2007h . Other authors have ap- 
plied conventional models for tidal damping to study 
the orb ital evolution of cl ose-in ex oplanets (Goldrcich 
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Ferraz-Mello et alj|2008l ; lJackson et al.ll2008l ). Although the 
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two processes may be co upled, few studie s have considered 
both processes together. iGu et al] (|2003T ) investigated the 
coupling between tidal heating and Roche lobe overflow, 
but they only considered the orbital effects of the tide 
raised on the planet by the host star. For close-in ex- 
oplanets on nearly circular orbits, like CoRoT-7 b, the 
tide raised on the host star by the planet can dominate 
the orbital evolution. In this section, we describe our 
model for coupled mass-orbital evolution and how we can 
use it to constrain the original orbit and mass of CoRoT-7 b. 
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2.1 Evaporative Mass Loss Model 



For energy-limited evaporative mass loss, the rate at which 
escaping gas molecules carry away energy from the planet is 
roughly propo rtional to the rate of input of energy from stel- 
lar insolation IjYelle et alj|20081 h Relating the rate of energy 
input from insolation to the change in gravitational energy 
required for a gas molecule to escape yields esti mates of mass 
loss r ates. We take the mass loss rate to be (jErkaev et al.l 
l200Tt ): 



dM p 
dt 



7Vrip€-T xuv 
GMpKtide 



(i) 



where F xuv is the extreme UV (XUV) flux from the star 
(wavelengths from 0.1 to 100 nm), evaluated at the planet's 
orbital distance. R p is the planet's radius, and G is the grav- 
itational constant. This equation assumes that the planet's 
optical cross- s ection in t he XUV is "k RJ, as suggested by 
lErkaev et all (|2007l ) and lYelld (|2004T ). The factors e and 
Kude (each of which may range from to 1) are discussed 
next. 

The parameter e represents the fraction of the incom- 
ing energy that is carried away by the escaping gas. Some of 
the incoming energy may drive chemistry or be lost through 
radiative cooling of the planet, rather than driving es cape. 
In these cases, e is less tha n 1 IjHubbard et a l. 2007?.;|YelleJ 
12004 ; iLammer et al.l 120091 ). For example, estimates of the 
mass loss rate for HP 209 458 b suggest —^f- ~4x 10 10 
g/s (e.g., lYelle et al.l 120081 ) . which corresponds to £ ~ 0.4. 
On the other hand, an analysis of atmospheric escape from 
Venus suggests 0. 1 5 may be a more appr opriate value for e 
|Chassefierel Il999 ) . ILammer et al.l (|2009j ) indicate that e is 
probably less than 0.6 for hot Jupiters but will depend on 
the host star's activity and the planet's atmospheric compo- 
sition, among other things. Moreover, vaporization and re- 
moval of Mercury's exosphere involve processes that may be 
differ ent from evaporation of gas giant planets l|Killen et al] 
120071 ). and evolution of the evaporating p lanet's composition 
may result in a time-dependent e- value (|Murrav-Clav et al] 
120091 ). However, Equation \T\ provides a reasonable estimate 
of evaporation of planetary mass by XUV flux, which plays 
an important, if not dominant, role for Mercury as well as 
gaseous planets. Here, we consider a range of time-invariant 
e values: 0, 0.1, 0.25, 0.5, and 1. We include e = to compare 
orbital evolution rates with and without mass loss. Note that 
we do not include any possible dependence of e on the solid 
planet's composition, even though mass loss may be more 
rapid if the planet were more enriched in volatiles. The de- 
tailed modeling required to quantify that effect is beyond 
the scope of this study. 

For planets close to their host stars, gas molecules may 
escape from a planet even with a velocity significantly less 
than the usual escape velocity because they only have to 
reach the Roche lobe. This reduction in the required veloc- 
ity results in an increased escape rate for a given rate of 
energy input. In Equation [T] Kude represents this reduction 
in required escape energy (jErkaev et al. . 2007). Kude is given 
by: 

3 



where £ is the ratio of the Roche radius (= (M p /3M*) 1//;i a) 
to R p . For a planet far from its host star, where the Roche 
radius is much greater than R p , Kude approaches one, indi- 
cating no enhancement of the escape rate. For HD 209458 
b, with a — 0.047 AU, Kude — 0.65. As tides bring a planet 
closer and closer to its host star, Kude decreases to reflect 
the reduction in required escape velocity. 

In employing this formulation for the mass loss rate, we 
treat its effects as if the mass loss is isotropic, that is, as if 
the flux of escaping gas is the same in all directions. In real- 
ity, gas may escape preferentially through the LI Lagrange 
point between the pl anet and star, where the net accelera- 
tion is close to zero (|Gu et alj 120031 ). which may affect the 
planet's orbital evolution. However, incorporating this di- 
rectionality may require a 3-D hydrodynamic model and so 
is beyond the scope of our study. We discuss this point more 
in Section 4. 

We also include the evolution of the stellar XUV flux 
in our model. The stellar XUV flux powers mass loss, but 
as stars age, typically their rotation rates drop with time, 
whic h reduces t heir magnetic activi ty an d XUV radiation 
{e.g., Skumanich 1972). Ribas et al. (2005) studied the evo- 
lution of the XUV flux from G-type stars like CoRoT-7 and 
suggested F xuv evolves according to: 



/ 

~ " u Gyr/ a 



= XtZ 
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where A = 29.7 ergs s _1 cm' 2 , f3 — 1.23, tG yr is the 
star's age in billions of years and a is in AU. This equa- 
tion is an empirical fit to data for observations of stars with 
ages greater than 0.1 Gyr and therefore may only apply for 
t > 0.1 Gyr. However, planets in our solar system (and, pre- 
sumably, in others) developed atmospheres in less than a 
few t ens of millions of years after formation l|Abe fc Matsuil 
Il986h . Also, before t = 0.1 Gyr, CoRoT-7's XUV flux was 
probably somewhat larger than at later times. Thus, atmo- 
spheric loss may also have been significant at earlier times 
than considered here. In this sense, our results provide a 
lower limit to the total mass lost by CoRoT-7 b. However, 
that earlier period of time was much shorter than the time 
span we consider, and so the total mass lost may have been 
not much greater than suggested by our calculations. With- 
out more data on stellar XUV at early times, stopping our 
models at t — 0.1 Gyr is probably an adequate approxima- 
tion for our purposes. 

CoRoT-7 b receives sufficient insolation that vaporiza- 
tion o f its rocky surface may p r oduce a tenuous atmos phere 
{e.g., ISchaefer fc Feglevl 120091 ) . IValencia et ail {|2009D sug- 
gested that the rate of vaporization of rocky material may 
be orders of magnitude larger than the atmospheric mass 
loss rate, sufficiently rapid to continually resupply the at- 
mosphere, so in our model, we allow evaporation to continue 
even after a model planet loses its entire original H/He en- 
velope. 

To determine the planet's radius after the original H/He 
envelop e is lost, we use Equations 7 & 8 from Fortnev et al 
2007b), along with corrections given in iFortnev et al 
2007aj). These equations provide mass-radius relations for 
solid planets made of ice and rock or rock and iron, which 
requires us to make some assumption about the planet's 
composition in order to determine its radius, given a mass. 

Figure Q] shows the range of compositions allowed by 
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CoRoT-7 b's observed mass and radius based on the equa- 
tions from lFortnev et al.l (|2007bft . That figure shows that, if 
made of ice and rock, CoRoT-7 b's ice mass fraction (IMF) 
may lie between 0.01 and 0.212. If made of rock and iron, 
CoRoT-7 b's rock mass fraction (RMF) may lie in the range 
0.719 to 0.983. These latter numb ers are similar t o comp osi- 
tions for CoRoT-7 b proposed by I Valencia et al,l i|2009l ). In 
our models, we consider the allowed minimum, middle and 
maximum values for both the IMF and RMF. However, the 
results are much less sensitive to our choice of IMF or RMF 
than to the other model parameters. 

Given the uncertainty in current composition and struc- 
ture of CoRoT-7 b, our model of coupled orbit evolution and 
evaporative mass loss cannot uniquely determine the mass 
of the gaseous envelope that CoRoT-7 b may originally have 
had. CoRoT-7 b may have begun as a gas giant with a core 
mass slightly larger than CoRoT-7 b's current mass (say, 
7 MEarth)- After evaporation stripped the original atmo- 
sphere, it could also remove ~ 1 MEarth from the original 
core to give the mass seen today. Consequently, in order to 
test the hypothesis that CoRoT-7 b started out as a gas gi- 
ant, we consider a range of masses for the core of the original 
planet (M core ), ranging from CoRoT-7 b's current mass (in 
which case, the planet very recently lost its original atmo- 
sphere) up to 19 MEarth- In our simulations, we run time 
backward, allowing the solid mass of the planet to increase 
until M p reaches our assumed M core . From that point con- 
tinuing back in time, we assume the planet "gains" mass in 
the form of H/He gas. 

Once the planet has a thick enough atmosphere, we de- 
fine the planet's radius R p to be the alt itude level at which 
most of the stellar XUV is ab sorbed (JErkaev et al.l 120071 ; 
lYellell2004l ; lLammer et al.ll2009ft . For M p ^ 0.018Mj„ p , R p 
corresponds to the 1 mbar pressure level. To determine R p 
as a function of M p and age for these plan ets, we use an 
update d version of the model described in iFortnev et al.l 
l|2007bl ). (This model provides R p at a series of grid points 
in M p , a and age, and we linearly interpolate between model 
grid points.) Although recent studies of evaporation of ex- 
oplanet atmospheres suggest tak ing R v at the 1 mbar level 
is the most appropriate choice (JYelld 12004 ). previous stud- 
ies suggest R v shoul d be taken at about the 0.1 nbar level 
l|Watson et al.l Il98ll ). This choice would significantly in- 
crease our calculated mass loss rate because R p would be 
much larger than what we've chosen, and our estimates here 
would provide lower limits on the total mass lost. 

For M p < 0.018 Mj up , the previous model for R p breaks 
down, so, for planets with such masses, we determine R v 
by assuming the planet's density remains constant at its 
value when M p — 0.018 Mj up . This assumption undoubt- 
edly introduces some error into the value for R p , but the 
relationship of radius to mass remains uncertain for planets 
transitioning from gas giants to rocky planets with atmo- 
spheres of several Earth masses. The atmospheric structure 
for these pla nets is a complicated function of atmospheric 
comp osition (|Miller-Ricci fc Fortnevll2010l : lRogers fe Seaeerl 
|2009j). Our approximation gives roughly the correct radii for 
Neptune and N eptune- like extra-so lar planets and has been 
applied before II Jackson et al.ll2008l ). Thus, it should suffice 
for the current study. We've also assumed that the mass 
loss has a negligible effect on the pla net's thermal evolut ion, 
consistent with previous work (e.o.. lTrilling et alJll99a ). 



2.2 Orbital Evolution Model 

To model the effects of tides on the orbit, we 
use a standard tidal model that has been applied 
in many studies of solar and extra-solar systems 



llGoldreichlll963l;lGoldreich fc Soterlll966l;lRasio et alj[l99^ ; 
iFerraz-Mello et all l20Qg| ; I Jackson et al.l 12008ft . Like many 
other close-in exoplanets, CoRoT-7 b raises significant tidal 
bulges on its host star. For planets that revolve more quickly 
on circular orbits than their host stars rotate, like CoRoT-7 
b, the tidal bulge on the host star lags behind the planet. 
This interaction results in a torque that spins up the host 
star and reduces the orbital semi-major axis a as: 

§ = -9 (G /M») 1/2 ^a-"/ 2 (4) 

where M* is the stellar mass, R* is the stellar ra- 
dius, and Q« is the modified tidal dissipation parameter 
for the star JGoldreich fc Soterl Il966l ; IJackson et all I2008J : 
lOgilvie fc Lin 12007). Because constraints on Q' t are poor 
{e.g., IJackson et al.l 120081 ), we consider a range of values: 
10 5 , 10 6 , and 10 7 . 

In using Equation [4] we assume CoRoT-7 b's orbital 
eccentricity has been negligible since t = 0.1 Gyr. CoRoT-7 
c might play a role in planet b's orbital evo lution by pump- 
ing up its orbital eccentricity ( [Barnes et alj|2010f ). However, 
tides probably keep CoRoT-7 b's eccentricity damped to 
small values. We also assume planet c has experienced neg- 
ligible mass loss and tidal evolution over its lifetime, due to 
its greater distance from the host star. Orbital evolution at 
that distance for a planet with a mass as sm all as C0R0T- 
7 c's minimum mass is also much smaller IjJackson et al.l 
l200Sft . 

To constrain CoRoT-7 b's initial mass (M Pt i n u) and 
semi-major axis (ainit), we numerically integrate Equations 
[T] and U together, backward and forward in time, and con- 
tinually evaluate all the other equations along the way. We 
consider the full range of current system parameters allowed 
by observational uncertainty, and run a unique integration 
for each combination of system parameters. That is, we run 
models with CoRoT-7 b's current mass M PtCur = 4.0, 4.8 
and 5.6 MEarth and with its current semi-major axis a cur — 
0.01691, 0.0172, and 0.01749 AU. Co RoT-7's age is r eport ed 
to lie between 1.2 Gyr and 2.3 Gyr (JQueloz et al.l l|2009h ). 
so we also run models with the planet's current age equal to 
1.2 and 2.3 Gyr. Along with the range of model parameters 
we also consider, altogether we run more than 10,000 unique 
calculations. 



3 RESULTS 

Our results show that the coupling of mass loss and orbital 
evolution may have played a significant role in CoRoT-7 b's 
history: Loss of planetary mass can reduce the rate of or- 
bital evolution because the rate of orbital evolution scales 
as the mass (Equation 2} . At the same time, even modest 
orbital evolution can significantly modify the mass loss his- 
tory of the planet because the mass loss rate depends so 
sensitively on the planet-star distance (Equations [TJ and [5]). 
In this section, we first consider the orbital and mass evo- 
lution if CoRoT-7 b has always been a solid planet. Then 
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we consider the evolution if CoRoT-7 b started out as a gas 
giant . 

3.1 Mass loss and orbital evolution if CoRoT-7 b 
has always been a rocky planet 

First we consider mass loss with negligible orbital evolu- 
tion. Figure [2] shows the mass evolution of CoRoT-7 b back- 
ward and forward in time, assuming a cur — 0.0172 AU, 
Mp lCU T = 5.6 M Ea rth, and R p = 1.6 R E arth, which corre- 
sponds to a constant rock mass fraction (RMF) = 0.719. 
At t — 0.1 Gyr, we see that CoRoT-7 b may have started 
with considerably more than its current mass. For exam- 
ple, if e = 1, M P: i„it could have been as large as about 9.2 
MEarth- Even for e as small as 0.1, M Pl mu was about 0.5 
MEarth larger than M P:Cur . As time moves forward, dM p /dt 
drops, largely as a result of the decrease in F xuv , and, by 
assumption, all mass evolution lines converge on M PtCur = 
5.6 MEarth at t — 2.3 Gyr. Though the current mass loss 
rate is smaller than in the past, for e = 1, dM p /dt ~7x 10 7 
kg/s (about 0.4 MEarth /Gyr). Obviously, smaller e- values 
give smaller mass loss rates. As time continues into the fu- 
ture, the mass loss rate continues to drop, and for e — 1, the 
planet will lose about 1 MEarth in a few billion years. 

In such proximity to its host star, tides likely do drive 
significant orbital migration. Figure [3] illustrates the mass 
evolution for the same a cur , M PtCUT , and (RMF) values as 
in the previous figure, but this time the calculation includes 
orbital evolution for a range of stellar Q' t . We can see from 
this figure that orbital migration plays an important role in 
both the planet's past and future. 

Starting again at t = 0.1 Gyr and moving forward in 
time, consider first the red lines, for which Q' t — 10 7 . In 
this case, CoRoT-7 b's orbit changes little, and ai n u is only 
about 0.001 AU larger than a cur . Thus, the planet receives 
more XUV radiation at all points in time than if Q' t were 
smaller. As a consequence, M Pt i n it is larger for larger Q'„- 
For example, for e — 1 and Q'„ = 10 7 , the planet loses about 
4 MEarth between t — 0.1 and 2.3 Gyr. By comparison, for 
e = 1 and Q' t = 10 , the total mass lost in the same period 
is about 3 MEarth- Thus, the total mass lost by CoRoT-7 b 
depends on Q' t and e. 

Now consider orbital evolution proceeding forward in 
time starting at t — 2.3 Gyr. Again, by assumption, the 
mass and orbital evolution lines pass through M P]Cur and 
a cur at t = 2.3 Gyr. CoRoT-7 b's orbit quickly decays, and 
the planet encounters the stellar surface and is destroyed 
in less than 2 billion years for Q' t ^ 10 6 . In this case, the 
planet only loses about 0.5 MEarth before it is destroyed. For 
Q' t — 10 7 , orbital evolution is much slower, and the planet 
loses more mass: for e = 1, about 1.25 MEarth over the next 
few billion years. Comparing the e = 1 lines in Figures[2]and 
[3] we see that the planet loses about 0.25 MEarth more by 
t = 10 Gyr with orbital evolution than without. Also, we see 
that different choices of e-values result in slightly different 
orbital evolution rates, as indicated by the divergence of the 
different linestyles in panel (b) of Figure [3] 

Given the short time CoRoT-7 b has left for Q'„ < 10 6 , 
we might conclude that Q'„ is much larger. Otherwise, it 
might seem unlik ely that we would see CoRoT-7 b at all. 
iHebb et all (|2010T l propose a similar argument for the planet 
WASP-19 b, which also seems on the verge of plummeting 



into its host star for Q' t ^ 10 9 . On the other hand, close-in 
exopla nets may well be com monly destroyed through orbital 
decay (| Jackson et al.ll2009f ). so CoRoT-7 b and WASP-19 b 
may just be next in a long line of planets to be destroyed. 

Although the amount of mass lost by CoRoT-7 b de- 
pends on Q' t and e, the planet's current orbit, mass and 
radius, as well as its composition, are also important for de- 
termining the total mass lost by the planet over its lifetime. 
However, by considering the full range allowed for these pa- 
rameters, we can constrain the change in mass and the ac- 
companying change in semi-major axis. 

Figure fj] shows the range allowed for the mass lost by 
CoRoT-7 b (AM P = M Pl i n it — M Pt cur) and the range allowed 
for ai n u, as functions of Q' t and e. In panel (a), the solid con- 
tours show the minimum value for AA/ P as a function of Q'„ 
and e, while the dashed contours show the maximum value. 
For exam ple, consider * = 10 and e = 0.4, values sug- 
gested by Ijackson et al.l (j2008f ) and the mass loss rate for 
HD 209458 b (Section 1), respectively. That point lies to 
the left of the "2 MEarth' dashed contour, which indicates 
the total mass loss cannot exceed 2 M Earth- The point also 
lies to the right of the "0.5 MEarth" solid contour, which 
indicates the total mass loss must exceed 0.5 MEarth- For 
the same Q' f and e-values, panel (b) tells us CoRoT-7 b's 
(Unit lay between 0.022 and 0.024 AU, at least 30 % larger 

Til till Q/cur • 

In interpreting Figure [3] panel (a), it's important to 
keep in mind that a dashed contour provides constraints on 
AM P only for the region to the left of the contour, while 
a solid contour provides constraints only for the region to 
the right of the contour. For example, consider the region 
sandwiched between the dashed "1 MEarth" and the solid 
"0.5 MEarth" contours. In this region, we can only conclude 
the change in mass in this region is less than 2 MEarth be- 
cause it lies to the left of the dashed "2 MEarth" contour. 
The contours nearest this region do not directly provide con- 
straints on this region. Otherwise, we might conclude para- 
doxically that, in this region, the total mass lost is greater 
than 1 MEarth and less than 0.5 MEarth- Similar consid- 
erations apply to panel (b), except that constraints from a 
dashed contour apply to the region below the contour, and 
constraints from a solid contour to the region above. 

Figure[4]also illustrates how the total mass lost depends 
on the orbital evolution and how, to a lesser extent, the or- 
bital evolution depends on the mass loss. In panel (a), for 
smaller Q' t , the contour lines become more horizontal. This 
trend reflects the fact that, for faster orbital evolution, the 
total mass lost by CoRoT-7 b depends more sensitively on 
the chosen Q'„-value, as well as on e. A similar but less pro- 
nounced trend is evident in panel (b). The slight upward 
slope of the lines reflects the fact that, for a fixed Q'„-value, 
larger e-values result in more orbital migration during the 
planet's history. 

The prese nce of CoRoT-7 c in a nearby orbit (with 
a = 0.046 AU lOueloz et all [2009]) would seem to provide 
some dynamical constraints on the allowed M Pt i n it and a,init- 
values for planet b. We expect that, whatever the original 
mass and orbit of CoRoT-7 b, gravitational interactions be- 
tween the planets should never have ejected either planet 
from the system or changed the ordering of their orbits, i.e., 
the system has always been dynamically stable. Unfortu- 
nately no simple criteria have been found that determines 
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whether a given sy stem is dy n amica ll y stable in this way. 
Howev er, based on iGladmanl (1 19931 ), iBarnes fc Greenberd 
(2006) showed that the boundary for this kind of stability 
for a two planet system can be approximated by the follow- 
ing relationship: 



( i V®\ , . «2 . , 04/3/^1/^2 

(/*i + -^) (Mi + HiS) >l + 3' ^^ 



(•») 



where fj,i is the ratio of planet i's mass to the total system 
mass, a — fix + jj,2, and 5 = (02/01 ) 2 . (We've assumed the 
orbits are circular.) Index 1 refers to the inner planet (b), 2 
to the outer planet (c). For planet masses and semi- major 
axes that satisfy the inequality, the corresponding system is 
likely to be dynamically stable. For planet masses and semi- 
major axes that do not satisfy the inequality, the stability 
is undetermined. Inequality (JSJ also does not apply to two 
planet systems in a low-order mean motion resonance. 

All M P A n it and a; n ; t -values represented in Figure [3] 
satisfy the inequality, which suggests the planetary system 
would have been stable for all cases indicated in the figure. 
However, our calculations do suggest that CoRoT-7 b may 
have passed through a 3:1 mean motion resonance, as its a- 
value passed through 0.022 AU. Inequality ([5} does not pro- 
vide any information about the stability of such a resonance 
encounter, so N-body modeling that includes tidal damping 
would be required to determine the stability. However, pre- 
vious studies of mean motion re sonances in satellite systems 
{e.g., (|Hamilton fc Zhandl2007l )) suggest that, if CoRoT-7 
b encountered such a resonance, the planet could migrate 
through it without being trapped (although the eccentricity 
could be temporarily boosted). Future studies should ad- 
dress such scenarios. 



3.2 Mass loss and orbital evolution if CoRoT-7 b 
is the remnant of a gas giant 

If CoRoT-7 b began instead as a gas giant planet, the to- 
tal mass lost may have been much larger. Recall, however, 
that determining the mass of the original atmosphere re- 
quires knowledge of the original core mass. Because, in our 
model, we allow the solid mass of the planet to evaporate 
after the original atmosphere evaporates, we can't estimate 
uniquely both the original core mass and the mass of the at- 
mosphere. To circumvent this problem, we consider a range 
of values for the original mass of the core M core . In the 
calculations described below, as we run time backward, we 
allow the solid mass of the planet to increase until it reaches 
the chosen M core . Going back in time past that point, we 
allow the atmospheric mass to increase until t — 0.1 Gyr. 
In our models, while M v ^ M core , we assume the planet is 
a gas giant, with an H/He envelope. Once M v drops below 
Mcore, we assume the planet has completely shed its H/He 
envelope and is only a denuded rocky/icy body. The critical 
core mass above which a rocky/icy body will begin to ac- 
crete a significant H/He envelope from the protoplanetary 
disk may span the ran ge from 1 to more than 10 MEarth 
fe.q.. llkomaet al.|[200il ). 

To illustrate the importance of our choice of M core , Fig- 
ure[5]shows mass evolution for two different assumed M core - 
values: 6 (red lines) and 19 MEarth (blue lines). The calcu- 
lations use the same model parameters as in Figures [2] and 
[3] and neglects orbital evolution. From t = 2.3 Gyr forward, 



the mass evolution for both lines is the same as in previous 
figures. In the past, though, there are important differences. 

Consider first the blue lines, for which we've set M core = 
19 MEarth- By assumption, while M p < M core , the planet 
is only a rocky/icy body with no H/He envelope. For this 
calculation, we see that the planet's original mass M Pi i n it 
was always less than 19MEarth- Thus, if we require that 
the planet's initial mass must have been at least 19MEarth 
to accrete a gaseous envelope from the protoplanetary disk, 
then this model planet never accreted such an envelope. As a 
consequence, the planet's radius R p is always much smaller 
than it would be if it managed to accrete an atmosphere. By 
comparison, for the M core = 6 MEarth case (red lines), the 
solid component of the planet was at least 6 MEarth about 
two billion years ago for all e ^ 0.1. Thus, according to our 
assumptions, the planet may have started with a significant 
gaseous envelope and, consequently, its radius would have 
been much larger during those times than if it never had 
an atmosphere, as for the previous model planet. As an ex- 
ample, for e — 1, the atmospheric mass could have been as 
large as about 200 M Earth- 

For the models we consider in this subsection, our choice 
of e plays a role that it did not play in the previous subsec- 
tion in determining the qualitative nature of past mass loss. 
As expected, for smaller e-values, the original mass of the 
planet was smaller. Also, because we assume the planet shed 
its original atmosphere before t = 2.3 Gyr, the smaller mass 
loss rates mean that the loss of the original atmosphere must 
have occurred further in the past. For example, for e = 0.5 
and 1.0, the model planet not only started with much larger 
atmospheric masses than for smaller e, the gaseous envelope 
was shed as recently as 500 million years ago. By contrast, 
for e = 0.25 and 0.1, the original atmosphere's mass was less 
than several MEarth and must have been shed at least 1.5 
billion years ago. 

However, we see that the lines for e = 0.25 and 0.1 stop 
before reaching t — 0.1 Gyr. The lines stop because, for 
slightly more massive atmospheres, the radii of these plan- 
ets at these early times would have exceeded the Roche lobe, 
and our mass loss model would break down. By contrast, 
the model planets with e = 0.5 and 1 are massive enough 
at these early times that their atmospheres remain bound. 
This dichotomy is a consequence of the thermal evolution 
of the planets' interiors. At early times, the interiors of our 
model gas giants are much warmer than at later times, and 
the planet s experience s i gnificant cooling and c ontraction as 
they age ( JGuillotl I200BI ; iFortnev et al.ll2007bh . Before t ~ 
1 Gyr, an atmosphere for a planet of only several MEarth 
would be so warm and the planet's surface gravity so small 
that the planet's radius would exceed the Roche limit. A 
much more massive planet has sufficient gravity that most 
of its atmosphere remains bound. 

What does all this mean for the history of CoRoT-7 b? 
For small mass loss efficiencies (small e) and negligible or- 
bital migration, CoRoT-7 b wouldn't have had enough time 
since its formation to lose more than a few Earth masses 
of gas through evaporative mass loss alone. Instead, if the 
planet started as a gas giant with a mass similar to Nep- 
tune's and orbital migration has been negligible, CoRoT-7 
b probably underwent a period of Roche lobe overflow early 
in its life. 

If CoRoT-7 b has undergone significant orbital migra- 
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tion, though, it may have avoided this early phase of Roche 
lobe overflow because it would have been farther from its 
host star than it is now. Figure [6] illustrates the planet's 
mass evolution including orbital migration for an assumed 
M core — QM Earth- Going into the future, we see results simi- 
lar to previous figures. There are again important differences 
in the past, though. 

As in Figure [3] moving forward in time from t = 0.1 
Gyr, orbital decay increases the insolation received by the 
planet, promoting mass loss. Consequently, for a fixed e- 
value, smaller Q'.-values correspond to smaller M Pt mu in 
general. For example, for Q' t = 10 (blue lines) and e ^ 0.25, 
mass loss rates are always small enough that the planet 
would not have had time to shed a large gaseous envelope 
by t = 2.3 Gyr. 

Figure [6] also shows the important feedback on the or- 
bital migration resulting from substantial mass loss. Panel 
(b) shows that, for the planets that begin with substan- 
tial gaseous envelopes, orbital decay abruptly slows once 
the planets lose their massive atmospheres. These results 
show that, if CoRoT-7 b has undergone enough orbital mi- 
gration, it may have avoided an early phase of Roche lobe 
overflow. The results also suggest that, in order to retain 
their atmospheres, any very close-in Neptune-like planets 
found orbiting sun-like stars must have started farther away 
from their stars. (Such restrictions may not apply to Gl 1214 
b because it orbits a muc h less massive star than CoRoT-7 
[Charbonneau et al.|[2009l ].') 

By running many such models, we can constrain the 
original mass and semi-major axis of a CoRoT-7 b that be- 
gan as a gas giant. Without knowledge of M core , though, 
we can only estimate the maximum value for M Pl i n it- Also, 
because Roche lobe overflow occurs for many of our model 
planets if they undergo little orbital migration, we can only 
meaningfully estimate the maximum value for ai n u as well. 

Figure [7] shows the results of such an analysis. (Because 
here we are only interested in the possibility that CoRoT-7 
b was a gas giant, we exclude models with e = for which 
there is no mass evolution.) In this figure, we see the same 
trends indicated in previous figures, although they are more 
pronounced in this figure as a result of the larger masses 
involved. For the nominal values Q', = 10 5 ' 5 and e = 0.4, 
from panel (a), we see CoRoT-7 b could have started with 
M Pl init no larger than 110 Msarth- For these same model 
parameters, we see in panel (b) ai„u ^ 0.028 AU. Thus 
CoRoT-7 b may have undergone significant orbital migra- 
tion over its lifetime. Holding Q'„ fixed at 10 5 ' 5 , we see that 
larger e- values allow larger M Pt i n it, but the contours flatten. 
This flattening indicates that an increased mass evolution 
rate is accompanied by an increased orbital evolution rate, 
which helps to slow the mass evolution - the same feedback 
seen in Figure [4] Holding e fixed at 0.4, we see that larger 
Q'.-values likewise allow larger and larger M Pt i n u- 

We can again apply our simple stability criterion, In- 
equality (J5]), to determine whether any of the M Pt i n it and 
Oimt values in Figure [7] can be ruled out on dynamical 
grounds. All values shown, however, satisfy the inequality, 
suggesting we cannot use dynamical considerations to rule 
out the possibility that CoRoT-7 b was a gas giant. Also, as 
in Figure [4j passage through mean motion resonances may 
have occurred. 



4 DISCUSSION 

Our study suggests CoRoT-7 b may have begun with a much 
larger mass than it has today, perhaps even with a substan- 
tial gaseous component. A more precise determination of 
CoRoT-7 b's composition may allow us to corroborate or re- 
fute the suggestion that the planet began as a gas giant. In 
particular, the large mass loss rates indicated by our calcu- 
lations suggest CoRoT-7 b may currently have an extended 
exosphere of vaporized material, observations of which may 
provide clues to the planet's composition. Our results also 
suggest the coupling between orbital evolution and mass loss 
may have shaped the observed population of close-in exo- 
planets. Although our calculations have elucidated the roles 
of orbital evolution and mass loss in the origin of CoRoT-7 
b, further work will be needed to test our assumptions and 
consider other potentially important physical processes. 

For CoRoT-7 b and other close-in low-mass planets, it 
may be possible to distinguish remnant co res of gas giants 
based on estimates of their masses and radii (|Ravmond et al.l 
2008). Many previous studies have developed constraints 
on the compositions of exoplanets made of solid materials 



based on estimates of their masses and radii (Valenc ia et al 



20071 ; iFortnev et all l2007bl ; ISotin et all 120071 ; lAdams et al 



2008). These studies show that determining the composi- 



tions of solid planets of end-member compositions (e.g., 
planets made entirely of ice or iron) is relatively easy but 
determining the compositions of planets with mixed compo- 
sitions is more difficult. Moreover, it may be hard to con- 
clude that a solid exoplanet is indeed a remnant core because 
the compositions of the core of a gas giant may vary from 
planet t o planet, in respon se to the various formation con- 
ditions ijPollack et al.lll996l ). 

In any case, if a close-in solid planet were found to 
have a substantial volatile component, it is unlikely to have 
formed in situ and is probably th e remnant of a more m as- 
sive planet, perhaps a gas giant l|Ravmond et al J 12008 ). In 
fact, such considerations st rongly suggest the recentl y dis- 
covered planet GJ 1214 b ([Charbonneau et al.l 120091 ) may 
be the remnant of a much larger planet. Given the planet's 
apparent density, recent studies suggest the planet may have 
an H/He envelope which is on ly 0.05% of the planet's mass 
l|Miller-Ricci fe Fortnevl l2010i ) or a water layer that is al- 
most 5 0% of th e planet's mass (Ro gers fe Sea ger 2009). Also, 
ICharbonneau et al.l (|2009l ) calculate that the current mass 
loss rate for GJ 1214 b may be sufficient to remove its entire 
gaseous envelope in about 700 Myrs, suggesting the planet 
probably started with a la rger mass. W ith its very close- 
in orbit (a = 0.0179 AU [Charbonneau et alJl20Qg |). the 
coupling between mass loss and orbital evolution from tidal 
damping would have played an important role in GJ 1214 
b's history. (We do not include GJ 1214 b here, in part, be- 
cause our model for stellar XUV radiation is not well-suited 
for M-dwarf stars, like GJ 1214. 



Recently IValencia et al.l (J2009J) considered the con- 
straints on CoRoT-7 b provided its mass and radius. They 
suggested that the available constraints are consistent with 
a wide range of compositions and internal structures, from a 
planet made entirely of silicates to one with a massive water 
vapor atmosphere. They also considered the possibility that 
CoRoT-7 b is the remnant of an evaporated gas giant and 
showed that the planet may have lost many Earth masses 
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of gaseous material, consistent with our results. That study 
did not include tidal migration, which we've shown may be 
important. 

Other recent studies cast doubt on the possibility 
that CoRoT-7 b started out as a Jupiter-like gas giant. 
iMurrav-Clav et al] (2009) suggested that complete evapora- 
tion of a Jupiter-mass planet's atmos phere may requ i re un - 
realistically large mass loss efficiencies. lLammer et al] (2009) 
suggested that, even if CoRoT-7 b began with a significant 
gaseous envelope, its original mass was probably less than 
Neptune's for the same efficiencies we've used here. How- 
ever, we find that CoRoT-7 b may be the remnant core of 
a gas giant planet with a mass as large as Saturn's even for 
relatively small mass loss efficiencies. The reasons for this 
disagreement probably lie in differences in the model de- 
tails. 

For example, lLammer et al] 1120091 ) considered a model 
very similar to ours (although their model for the stellar flux 
was slightly different), but did not apply their model to a 
planet as close to its star as CoRoT-7 b, which receives more 
insolation and for which the "Roche lobe effect" is more im- 
portant. Thus it may ex perience more r apid m ass loss than 
the planets modeled by lLammer et al] |2009l ). In any case, 
it will be interesting to reconsider the results from previous 
studies by including the effects of orbital evolution as we 
have done. 

If vaporization of CoRoT-7 b continues today, CoRoT- 
7 b may b e surrounded by a dens e exosphere, similar to 
Mercury's. ISchaefer fc Feglevl (|2009f ) found thermal vapor- 
ization of rocky material could produce a column abundance 
of sodium of perhaps 10 19 crn~ 2 , almost 10 8 greater than th e 
abundance of sodium for Mercury (e.q.. lKillen et alj|2007h . 
For comparison, the column abundance of the Earth's atmo- 
sphere is ~ 10 25 cm~ 2 . Other processes that generate Mer- 
cury's exosphere (e.g., photon-stimulated desorption) would 
probab ly produce an eve n large r column abundance. Al- 
though ISchaefer fc Feglevl l|2009l ) did not include processes 
that might quickly remove vaporized materials from the ex- 
osphere (such as photoionization and stellar wind entrain- 
ment), such large column abundances for sodium and other 
materials suggest an exosphere around CoRoT-7 b might be 
observable. Even if such observations are not possible for 
CoRoT-7 b, they may be possible for other close-in rocky 
exoplanets soon to be discovered. 

Our results also suggest very close-in Neptune- like plan- 
ets orbiting Sun-like stars may only be able to retain their 
atmospheres if they have undergone significant orbital decay 
and started out farther from their host stars. As our under- 
standing of tidal dissipation within stars improves (see, e.g., 
IPenev et alj|2009l) . we may find that Q*-values significantly 
exceed the values considered here. Such a result might sug- 
gest we will be unlikely to find very close-in Neptune-like 
planets because those planets would have rapidly shed their 
atmospheres through Roche lobe overflow early in their his- 
tories. 

Our results have implications for the population 
of observed c l ose-in exoplanets as well. For example, 
iHubbard et al] l|2007r ) compared the mass distribution of 
close-in exoplanets to that of exoplanets far from their host 
stars and suggested they are statistically identical, and that 
the mass distribution expected if significant evaporative 
mass loss from close-in exoplanets had occurred is incon- 



sistent with the observed distribution. However, the sam- 
ple size of that study was limited (the number of close-in 
exoplanets has more than doubled since that paper was 
written), and the disagreement between observations and 
the expectations from mass loss models was marginal (less 
than t wice the observational sta ndard deviation). More re- 
centlv. lDavis fc Wheatlevl (|2009l ) noted a lack of planets with 
small masses and short orbital periods and suggested evap- 
oration has removed these planets. However, both of these 
previous stud i es neg lected the effects of orbital evolution. 
Ijackson et al] (|2009T ) noted a similar gap in the distribution 
of planetary semi-major axes and ages and found the gap is 
consistent with tidal destruction of planets, but that study 
ignored the evaporative mass loss that would accompany the 
orbital evolution. 

Our results here suggest that the most massive close- 
in exoplanets lose little mass, but will experience the most 
rapid orbital evolution. Smaller planets will lose more mass, 
but experience little orbital evolution. Together these pro- 
cesses may cause the most massive planets to quickly crash 
into their host stars without losing much mass. The least 
massive planets could be completely evaporated without 
much orbital evolution. Planets with intermediate masses 
may experience both significant orbital evolution and mass 
loss. Such a scenario would undoubtedly influence the mass 
distribution for close-in exoplanets. However, the exact bal- 
ance between loss and replenishment could be complicated 
and merits further study. 

Our results have shown the importance of considering 
tides when studying mass loss of close-in planets and partic- 
ularly the role tidal considerations can play in constraining 
the original masses of close-in planets. However, a number 
of issues remain to be addressed. For example, the depen- 
dence of R p on mass for planets with gaseous envelopes of 
only several Earth masses is uncertain. For these planets, 
R p likely depends in a complex way on the envelope's com- 
position, but the composition of such planets is poorly con- 
strained and may span a wide range. In fact, for such planets, 
of which GJ 1214 b may be an example, knowledge of the 
planet's dens ity is not sufficient for determining the planet's 
comp osition (|Rogers fc Seager 2009; M iller-Ricci fc Fortnevl 
120101 ). Moreover, significant atmospheric mass loss from gas 
giants likely results in a time- variable atmospheric composi- 
tion, as the lightest constituents may be lost first. Studying 
these details may be especially important for understanding 
the origins of exoplanets with low masses, which will likely 
be found in close-in orbits. 

Another potentially important effect we've neglected 
is the channeling of escaping gas by the stellar gravity. 
The tidal gravity of the host star may cause mass to 
flow preferentially throu gh the LI point l|Gu et al] 120031 ; 
IMurrav-Clav et al] [2009) . This effect would be especially 
important during Roche lobe overflow, when a substantial 
portion of the planet's atmosphere may become unbound 
from the planet. Exchange of momentum between the escap- 
ing gas and the remaining planet can counteract the inward 
pull of tides, even c ausing the planet to recede from the star 
|Chang et alj|2010l ). This effect may significantly modify the 
orbital evolution we've considered here and may be responsi- 
ble, at least in part, for the clustering of th e orbital periods 
of clo se-in exoplanets around three days (|Cumming et al] 
|200S| ). 
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5 CONCLUSIONS 

The first discovered rocky exoplanet, CoRoT-7 b, could 
actually be the remnant core of a gas giant. The planet's 
proximity to its host star and strong stellar insolation 
suggest that the planet may have lost a significant amount 
of mass over its lifetime. The proximity also suggests tides 
have played an important role in the planet's history, 
causing it to migrate inward. By coupling models of 
evaporative mass loss and tidal evolution, we find that the 
interplay between orbital decay and mass loss can lead to 
a wide variety of evolutionary pathways. Using a model 
that couples these two processes, we obtain constraints on 
the masses and orbits of gas giants that may have given 
rise to the CoRoT-7 b we see today. Our results suggest 
a complex orbital history for the planets in the CoRoT-7 
system, possibly including passage through mean-motion 
resonances, even if CoRoT-7 b has always been a rocky 
planet. 

Our calculations have general implications for the his- 
tories of both close-in rocky and gaseous planets. Previous 
studies have pointed out the signatures of evaporation and 
tidal evolution on the distribution of masses and semi- major 
axes for gaseous close-in exoplanets, but these studies have 
not considered the coupling of the two processes. Although 
our calculations allow CoRoT-7 b to be the remnant rocky 
core of a gas planet, confirming this hypothesis will require 
additional evidence. For example, improved constraints 
on the compositions of gas giant cores might allow us 
to distinguish remnant rocky cores on the basis of their 
physical properties. More detailed dynamical studies of 
CoRoT-7 and other similar systems will provide important 
constraints on their origins and histories. Whatever the 
details, mass loss and tidal evolution both probably played 
key roles in the histories and origins of many close-in 
exoplanets, rocky or gaseous. 
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Figure 1. Range of planetary compositions consistent with CoRoT-7 b's mass M p and radius R p , as reported in lLeger et al.l (2009). 
The filled circle represents the nominal M p and R p , while the shaded rectangle shows the uncertainties. Solid lines represent the M p -R p 
relation if the planet is made of rock and/or iron, while the dashed represent the relation if the planet is made of rock and/or ice. From 
bottom to top, the solid lines represent rock mass fractions (RMF) of (pure iron), 0.719, and 0.983. From top to bottom, the dashed 
lines represent ice mass fractions (IMF) of 1, 0.212, and 0.0108. 
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Figure 2. Evolution of M p for CoRoT-7 b, neglecting orbital migration. The empty circle represents an assumed current mass M PtCur = 
5.6 M Earth and a current age of 2.3 Gyr. The different lincstyles correspond to different assumed evaporation efficiencies e, as indicated. 
For this figure, we've determined R p as a function of M p assuming the planet is composed only of rock and iron, with a rock mass 
fraction (RMF) = 0.719. 
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Figure 3. Evolution of M p and a for CoRoT-7 b, including orbital migration. The empty circle in panel (a) represents an assumed 



2.3 Gyr. Again, we've assumed (RMF) remains constant at 0.719. The different line colors correspond to different assumed values of Q'„ , 
as labeled. The solid black line represents the stellar surface. 
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Figure 4. Panel (a): Minimum and maximum values for AM P = M p i„n — M PlC ur found by considering the range of system and model 
parameters, as discussed in Section 2. Solid contours show minimum values for AM P , and dashed contours show maximum values. For 
example, AM P lies between 0.5 and 2 MEarth for e = 0.4 and Q, = 10 . Panel (b): Minimum and maximum values for ai„u, again 
by considering a range of parameters. Solid and dashed contours again represent minimum and maximum values, as in panel (a). For 
example, ai n u lies between 0.02 and 0.026 AU for e = 0.4 and Q'„ = 10 5 ' 5 . 
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Figure 5. Evolution of M p , neglecting orbital migration, for different assumed core masses M core , as indicated by the different line 
colors. The empty circle has the same meaning as in previous figures. While M p < M core and the planet is losing mass in the form of 
rocky material, we assumed a constant (RMF) = 0.719. 
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Figure 6. Evolution of M p and a for CoRoT-7 b, including orbital migration. Here, we've assumed (RMF) = 0.719 and M cor 

M Earth- 
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Figure 7. The maximum values allowed for CoRoT-7 b's initial mass M p j n jt and initial semi-major axis ai n i t . Dashed contours have 
meanings similar to the dashed contours in Figurel4l For example, M p i n u < HOJVfBorth an d a init < 0.03/1(7 for t = 0.4 and Q'„ = 10 5 ' 5 , 
but the contours provide no constraints on the minimum initial values. 



